Spectral and optical properties in the antiphase stripe phase of the cuprate 
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We investigate the superconducting order parameter, the spectral and optical properties in a 
stripe model with spin (charge) domain-derived scattering potential Vs (Vc). We show that the 
charge domain-derived scattering is less effective than the spin scattering on the suppression of 
superconductivity. For Vs 3> Vc, the spectral weight concentrates on the (7r,0) antinodal region, 
and a finite energy peak appears in the optical conductivity with the disappearance of the Drude 
peak. But for Vs ~ Vc, the spectral weight concentrates on the (7r/2,7r/2) nodal region, and 
a residual Drude peak exists in the optical conductivity without the finite energy peak. These 
results consistently account for the divergent observations in the ARPES and optical conductivity 
experiments in several high-Tc cuprates, and suggest that the "insulating" and "metallic" properties 
are intrinsic to the stripe state, depending on the relative strength of the spin and charge domain- 
derived scattering potentials. 
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I. INTRODUCTION 



The nature of spin and/or charge inhomogeneities, es- 
pecially in the form of stripes, in some cuprates and their 
involvement to high-temperature superconductivity are 
currently debate issuesii The stripe state is character- 
ized by the self-organization of the charges and spins 
in the Cu02 planes in a peculiar manner, where the 
doped holes are arranged in one-dimensional (ID) lines 
and form the so-called " charge stripe" separating the an- 
tiferromagnetic domains. The stripe-ordered state mini- 
mizes the energy of the hole-doped antiferromagnetic sys- 
tem, thus leading to an inhomogeneous state of matter. 
Static one-dimensional charge and spin stripe order have 
been observed experimentally in a few special cuprate 
compounds, specifically in Lai.6-£i;Ndo.4Sr2;Cu04^''^ and 
La2_2;Baa:Cu04 with x = 1/84*^ Similar signatures iden- 
tified in La2-a;Sr^Cu04 (LSCO)^^^!^ and other high 
temperature superconductorsiSiliii^ point to the possible 
existence of stripes, albeit of a dynamical or fluctuating 
nature. 

A pivotal issue about this new electronic state of mat- 
ter concerns whether it is compatible with superconduc- 
tivity, and possibly even essential for the high transition 
temperatures, or it competes with the pairing correla- 
tions. A prerequisite for addressing these issues is to un- 
derstand the electronic structures of various stripe states 
in different cuprates, and to answer the question whether 
the stripe phase is intrinsically " metallic" or " insulating" , 
given its spin- and charge-ordered nature. Angle-resolved 
photoemission spectroscopy (ARPES) study by Zhou 
et al. in (Lai.28Ndo.6Sro.i2)Cu04 with static stripes 
have found the depletion of the low-energy excitation 
near the (7r/2,7r/2) nodal region.'^ In another compound 
Lai.875Bao.i25Cu04, a system where the superconductiv- 
ity is heavily suppressed due to the development of the 
static spin and charge orders. Valla et al. have detected 



the high spectral intensity of the low-energy excitation in 
the vicinity of the (7r/2,7r/2) nodal region [while antin- 
odal low-energy quasiparticle near (0,7r) are gapped] 
The compound (Lai.4_a;Ndo.6Sra;)Cu04 (a:; = 0.10 and 
0.15) with static one-dimensional stripe, seems to be an 
in-between system, in where the existence of spectral 
weight around the nodal region, though weak, has been 
identified.^'* 

Meanwhile, optical conductivity measurements on 
the systems with a stripe phase also display the di- 
vergent results. In Lai.275Ndo.6Sro.i25Cu04*^ and 
Lai.875Bao.i25-£cSr2;Cu04fi2, a finite frequency absorp- 
tion peak with almost disappearance of the Drude mode 
in the low-frequency conductivity in several experiments 
has been interpreted as collective excitations of charge 
stripes or as charge localization from the disorder created 
by Nd or Ba substitutions. These observations may sup- 
port the suggestion that such stripe-ordered state should 
be "insulating" in nature.*^ On the other hand, opti- 
cal experiment on Lai.875Bao.i25Cu04*^ has observed 
a residual Drude peak with a loss of the low-energy 
spectral weight below the temperature corresponding to 
the onset of charge stripe order, which indicates that 
stripes are compatible with the so-called nodal-metal 
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Although, there have been some theoretical studies on 
the spectral and optical properties in the stripe phase 
in the past years ^^4, 25, 26, 27 ^j^g contradictory observa- 
tions in recent experiments as mentioned above have yet 
not been explained consistently in theoretical frame by 
adopting a realistic stripe model. In this paper, by using 
a stripe model in which the experimentally observed spin 
and charge structures at 1/8 doping are well reflected, 
we show that the spin domain-derived scattering will de- 
press the zero-energy spectral weight around the nodal 
regions, while the charge domain-derived scattering will 
suppress mostly those around the antinodal regions and 
the hot spots. Compared to the ARPES data, this sug- 
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gests that the different spectral weight distribution may 
result from the different relative strength of the spin 
and charge domain-derived scattering potentials inher- 
ently existing in these compounds. Meanwhile, a finite 
frequency peak in the optical conductivity appears with 
the disappearance of the Drudc peak in the case of the 
dominant spin domain-derived scattering. While, when 
the charge domain-derived scattering is comparable to 
the spin one, a residual Drude peak exists with the dis- 
appearance of the finite energy peak. This suggests that 
both the "insulating" and "metallic" properties are in- 
trinsic to the stripe state without introducing another 
distinct metallic phase. 

The rest of this paper is organized as follows. In Sec. 
II, we introduce the model Hamiltonian and carry out 
the analytical calculations. In Sec. Ill, we present the 
numerical calculations and discuss the results. In Sec. 
IV, we present the conclusion. 

II. THEORY AND METHOD 

As the above discussed compounds have a doping den- 
sity at or near 1/8, we will in this paper consider the 1/8 
doping antiphase vertical stripe state. A schematic illus- 
tration of its charge and spin pattern is presented in Fig. 
1. The charge stripes, with a unit cell of 8 lattice sites 



and 



Afc = 



As for the tight-binding energy band, we will choose the 
following ioxm^^^ 

ek= — 2{6t + J xo) {cos kx + cos ky) 

~46t cos kx cos ky — /i. (4) 



(Note for 1 /8 doping, there is one hole for every two sites 
along the length of a charge stripe) , act as antiphase do- 
main walls for the magnetic order, so that the magnetic 
unit cell is twice as long as that for charge order. Due to 
the periodical modulation of the stripe order, the elec- 
trons moving in the state will be scattered by the mod- 
ulation potentials. After Fourier transformation, the po- 
tential Vn can be written as the scattering term between 
the state k and those at k±nQ with Q = (37r/4,7r). Fol- 
lowing Ref. m, we expect that the terms Vi and V2 will be 
the dominant spin and charge domain-derived scattering 
term, and will be relabeled as Vg and Vc in the follow- 
ing, respectively. The weaker higher harmonic terms will 
be neglected here. In the coexistence with the super- 
conducting (SC) order, the model Hamiltonian can be 
written as a 16 x 16 matrix for k in the reduced Brillouin 
zone, 

iy^^'(7t(fc)(^^^^ ^^^^Cik), (1) 

where, the prime denotes the summation over the re- 
duced Brillouin zone. Ck is a column vector with its 
elements Ci{k) = Ck+(i~i)Q.^ for i = 1,2,- •• ,8, and 
Cl^,_^-_g^Q for i = 9, 10, • • • , 16. Both Hk and are 
8x8 matrix with 



(2) 



(3) 



I 

where, S is the doping density, and a d-wave SC order pa- 
rameter Ak — 2J Ao(cosfca; — cos ky) is assumed. Gener- 
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ally, the charge modulation will induce the modulation of 
the SC order leading to the finite momentum pairs. How- 
ever, in the present study, one of our aim is to examine 
the effect of the spin (charge) domain-derived scatter- 
ing on the SC order. In this regard, the average value 
of the SC order parameter is relevant and the modula- 
tion of the SC order will be ignored. We have checked 
the effect of this modulation and found no qualitative 
change in the results presented in Fig. 2. In the follow- 
ing, J = lOOmeV is taken as the energy unit, t = 2J, 
t — — 0.45t, J =1^- This dispersion can be de- 
rived from the slave-boson mean-field calculation of the 
t — t — J mode P^i^° , and in this way the parameters 
Aq, xo and /i are determined self-consistently. Here we 



Ao = 



Xo = 



n = 



where, Em{k) is the eigenvalue of the Hamiltonian, 
Umn{k) the elements of the matrix U{k), and f[E„i{k)] 
is the Fermi-Dirac distribution function. 

Then, the single particle Green functions (fc, iuJn) = 
— dr exp^'^"^'^ {TrCi{k, iT)Cj{k, 0)) can be expressed as 

and the spectral functions is 

A^{k,uj) = -ilmGij(fc,cj + ■iO+). (7) 



III. RESULTS AND DISCUSSION 

A. Self-consistent calculation of the SC order 
parameter 

We first present in Fig. 2 the self-consistent results of 
the SC order parameter as a function of and Vc. While 
the scattering from both spin and charge domain-derived 
scattering potentials in the stripe state leads to the sup- 
pression of the SC order parameter, the charge domains 
are more compatible with superconductivity than spin 
domains, as can be seen from Fig. 2(a). This may sup- 
port the statement that the SC pairing in the stripe state 



take it as a phenomenological form. In a self-consistent 
calculation, the Hamiltonian is first diagonalized by a 
unitary matrix f/(/c) with a set of trial values of Aq, xo 
and /i for given potentials Vg and Vc- Then Ao, Xo and 
/i are self-consistently calculated by using the relations: 
±Ao — (ci|C,;+T| — CiiCi+rt) (To get the d-wave pair- 
ing, the sign before Aq takes -I- for r = ztx and — for 
T = ±y, where x and y denote the unit vectors along 
X and y directions, respectively.), xo — X]cr('^L'^jo-), and 
n — X]cr(''L^*<T)' respectively. Reformularization of the 
expressions of Aq, xo and [i in terms of eigenfunctions 
and eigenvalues of the Hamiltonian, one obtains the self- 
consistency relations 



(5) 



occurs most strongly within the charge stripes, '^"'^ On the 
other hand, an interesting feature is that the SC order 
parameter will be zero at the spin domain-derived scat- 
tering potential Vsc ~ 0.14 in the absence of the charge 
domain-derived scattering, however, it will develop a no- 
ticeable value after turning on the charge domain-derived 
scattering potential, as shown in Fig. 2(b). This shows 
that the charge domain-derived scattering will lead to the 
emergency of the SC order which is otherwise destroyed 
by the spin only scattering. 

B. Distribution of spectral weight 

In Fig. 3, we present the distribution of the low-energy 
spectral weight in the original Brillouin zone (integrated 
over an energy window Ae — O.IJ about tp) in the 1/8 
antiphase stripe state for different spin (charge) domain- 
derived scattering potential Vg (K)- Let us first look 
at the limit where only the spin domain-derived scatter- 
ing is included, i.e., K = with = 0.15, one will 
find that the spectral weight around the nodal region 
is suppressed heavily[See Fig. 3(a)]. At another limit 
where only the charge domain-derived scattering is in- 
cluded {Vc — 0.17 with Vs = 0), the spectral weight 
around the nodal region is recovered and those around 
the hot spot (the cross of the Fermi surface with the line 
kx i ky = ±7r) and near the antinodal region are sup- 
pressed[See Fig. 3(b)]. Starting from the limit of 14 = 



J 



-'j;^^{cosk.^- cos ky) ^ Ui^{k)ulg{k)f[E^{k)] 

k m—1 
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and fixing Vs = 0.15, the spectral weight will redistribute 
gradually from the antinodal region to the nodal region 
with the increase of the charge domain-derived scatter- 
ing potential Vc, as shown in Figs. 3(c) and (d). When 
two scattering potentials are comparable, the strongest 
spectral weight situates around the nodal region, and at 
the meantime noticeable spectral weights along the whole 
Fermi surface is presented. Therefore, the divergent fea- 
tures observed in ARPES measurements by Zhou et al. 
in (Lai.28Ndo.6Sro.i5)Cu04^''- in which the low-energy ex- 
citations near the nodal region are depleted, and by Valla 
et al. in Lai.875Bao.i25Cu04i^ in which the high spec- 
tral intensity of the low-energy excitation in the vicinity 
of the nodal region is detected are consistently repro- 
duced here by a change of the relative strength between 
the charge and spin domain-derived scatterings. This 
consistent accounting enables us to propose that the spin 
domain-derived scattering dominates over the charge one 
in the former system while the scattering strengthes of 
them are comparable in the latter system. 

In the presence of the spin (charge) domain-derived 
potential, quasiparticles near the Fermi surface will be 
scattered from k to k ± nQ (n=l for the spin domain- 
derived potential, n=2 for the charge one), for the 1/8 
antiphase vertical stripe configuration shown as Fig. 1. 
This gives rise to two scattering channels from the spin 
domain with potential Vs, 



k 
k 



k + Q 
k-Q 



k 
k 



(3^/4, tt), 
(5^/4, tt). 



(8) 



and two scattering channels from the charge domain with 
potential Vc, 



k 
k 



k + 2Q = k 
k - 2Q = k 



(37r/2,0), 
(7r/2,0). 



(9) 



Strong potential scattering will destruct those parts of 
the Fermi surface connected by the above mentioned scat- 
tering wave vectors. Because the scattering wave vectors 
Q and —Q are close to the transferred momenta from 
the node to node scattering, so it will lead to a depletion 
of the spectral weight near the nodal region as shown in 
Fig. 3(a). On the other hand, the scattering wave vectors 
2Q and —2Q, which is near the connecting wave vectors 
between the two approximately parallel segments of the 
Fermi surface near the antinodal and hot spot region, 
the scatterings with these wave vectors will suppress the 
spectral weights around the antinodal and hot spot re- 
gions [Fig. 3(b)]. 



C. In-plane optical conductivity 

Now, we turn to the discussion of the in-plane op- 
tical properties in the 1/8 antiphase stripe state, and 
to see how they are influenced by the scattering from 
the spin and charge domains. We will fix the temper- 
ature at T = 0.05 in all calculations, in order to avoid 



the influence from the temperature induced change in 
the scattering rate. We consider an electric field ap- 
plied in the x direction, which is perpendicular to the 
stripe. From the Kubo formula for the optical con- 
ductivity, the real part of the optical conductivity is 
ai{Lu) = — limg_,o Iiii[n(q, w)]/w. The imaginary part of 
the current-current correlation function Im[n((j' —>■ 0,uj)] 
is given by 

16 

lm[n(g-0,c.)]= ^Y^'Y^yn^kyif^k) 

k j,i=i 

X / di^'ifiu + u)- fiio')] 



xAji{k,uj )Aij{k,LU + LU ). (10) 

Here, v^^ (k) is the diagonal element of the quasiparticle 
group velocity in the matrix form 







vik) = ( ^ 







(11) 



Figs. 4(a)-4(d) show the results for the optical conduc- 
tivity calculated with the same scattering potentials as 
used to get Fig. 3(a)-3(d). With only spin domain- 
derived scattering[Fig. 4(a)], no Drude-like component 
appears at zero frequency in the optical conductivity, in- 
stead a finite frequency conductivity peak occurs around 
0.3. This indicates that the system exhibits the "in- 
sulating" property— When only charge domain-derived 
scattering is considered [Fig. 4(b)], the Drude-like peak 
shows up and at the meantime the finite frequency peak 
remains. Optical conductivity involves the contribution 
from the quasiparticle excitations along the whole Fermi 
surface weighted by the quasiparticle group velocity. Due 
to the relative flat band structure near the antinodal re- 
gion for the high- Tc cuprates, the zero frequency optical 
conductivity mainly comes from the quasiparticle excita- 
tions around the nodal region. In the case of only spin 
domain-derived scattering, the nodal region of the Fermi 
surface is gapped and therefore the quasiparticle spectral 
weight is suppressed around the nodal region as shown 
in Fig. 3(a), so that the zero- frequency Drude-like peak 
is absent and a flnite frequency peak with its position 
being equal to the gap (w 2Vs = 0.3) occurs. For the 
charge domain-derived scattering, the gap opens around 
the hot spots and near the antinodal, but a large spec- 
tral weight situates around the nodal region, as can be 
seen from Fig. 3(b). Thus, the Drude-like peak emerges 
and the finite frequency peak remains (it is now situ- 
ates at « 2Vc = 0.34). As shown in Fig. 3(c), with 
the increase of the charge domain-derived scattering Vc, 
the gap near the nodal region which is resulted from the 
spin domain-derived scattering will be suppressed grad- 
ually and correspondingly the spectral weight will be en- 
hanced. As a result, the finite frequency peak in the opti- 
cal conductivity is shifted to lower frequency and the zero 
frequency component is lifted up gradually [Fig. 4(c)]. 
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When the charge domain-derived scattering is compara- 
ble to the spin one, the quasiparticles have noticeable 
spectra weight along the entire Fermi surface with its 
largest weight around the nodal region[Fig. 3(d)], then 
the Drude-like mode occurs at the zero frequency, and 
the finite frequency peak fades away and merges into the 
Drude-like peak, as shown in Fig. 4(d). The calculated 
results for the optical conductivity presented in Figs. 4(c) 
and 4(d) are consistent well with the experimental obser- 
vations in the stripe state of Lai.275Ndo.6Sro.i25Cu04ii^ 
and Lai.875Bao.i25Cu04i^, respectively. 



D. Discussion 

We now discuss the implication of our theoret- 
ical results. As noted in the introduction, in 
L^i 275Ndo gSro.i25Cu04 system, ARPES experiment 
has found that there is little or no low-energy spectral 
weight near the nodal region;^ and optical conductivity 
experiment has observed a finite frequency peak with al- 
most the disappearance of the Drude mode, indicating 
an "insulating" stripe state J^ii^ These spectroscopic fea- 
tures can be reproduced here with a strong spin domain- 
derived scattering potential Vs = 0.15 and a weak charge 
domain-derived potential Vc = 0.08 and 14 = 0, as shown 
in Figs. 3(a), 3(c), 4(a) and 4(c). Interestingly, in this 
parameter regime for the spin and charge domain-derived 
scattering, the SC order is destroyed as can be seen from 
Fig. 2(b). This is in consistent with the experimental 
fact that Lai.275Ndo.6Sro.i25Cu04 is nonsuperconduct- 
ing. In another cuprate Lai.875Bao.i25Cu04, ARPES 
spectra have identified the existence of high spectral in- 
tensity around the nodal regioUj-- and the optical con- 
ductivity measurement has observed a residual Drude 
peak without the finite frequency peakfi^ pointing to a 
so-called nodal metal state j-^^i-^^'^°i^-^'^^i^^ When compa- 
rable spin and charge domain-derived scattering poten- 
tials are assumed such as Vg — 0.15 and Vc = 0.17, we can 
reproduce these features consistently, as shown in Figs. 
3(d) and 4(d). On the other hand, a weak supercon- 
ductivity emerges in the otherwise nonsuperconducting 
regime (when only spin scattering potential Vsc is con- 
sidered) with the increase of the charge domain-derived 
scattering potential[see Fig. 2(b)]. This suggests that the 
weak superconductivity in Lai.875Bao.i25Cu04 is likely 
beneficial from the metallic behaviors of the stripe state 
originated from a sufficient charge domain-derived scat- 
tering. The above mentioned consistent accounting for 
both divergent spectroscopic features observed in two 
families of high- Tc cuprates indicates that the stripe 
state may be intrinsically " insulating" or " metallic" , de- 
pending on the relative strength of the spin and charge 
domain-derived scattering potentials. Specifically, a large 
spin domain-derived scattering potential favors the "in- 



sulating" state, while a large charge domain-derived scat- 
tering potential the "metallic" state. 

IV. CONCLUSION 

We have calculated the SC order parameter, the spec- 
tral function and the optical conductivity in a stripe 
model with spin and charge domain-derived scattering 
potentials {Vs and Vc). The self-consistent calculation of 
the SC order parameter shows that the charge domain- 
derived scattering is less effective than the spin scattering 
on the suppression of superconductivity, and may even 
lead to the emergency of the SC order which is otherwise 
destroyed by the spin only scattering. For Vs Vc, the 
zero-energy spectral weight disappears around the nodal 
points, and a finite energy peak appears in the optical 
conductivity with almost the disappearance of the Drude 
peak. But for Vs ^ Vc, the spectral weight concentrates 
on the nodal region, and a residual Drude peak exists in 
the optical conductivity without the finite energy peak. 
These results consistently account for the divergent spec- 
troscopic properties observed experimentally in two fam- 
ilies of high-Tc cuprates, and demonstrate that both the 
" insulating" and " metallic" behavior may be the intrinsic 
properties of the stripe state, depending on the relative 
strength of the spin and charge domain- derived scatter- 
ing potentials. 
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FIG. 1: (Color online) Schematic illustration of the charge 
and spin patterns in the 1/8 doped antiphase stripe state. 
Circles represent the charge domain wall (An empty circle 
indicates a hole density of one per site), and arrows the copper 
spins. 
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FIG. 2: (Color online) (a) Superconducting order parameter 
as a function of K and Vc, respectively, (b) A two-dimensional 
map of the superconducting order parameter in the parameter 
space of Vs and Vc- 
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FIG. 3: (Color online) Spectral weight distribution for differ- 
ent spin (charge) domain-derived scattering potentials in the 
normal state with (a) Vs = 0.15 and Vc — 0, (b) Vs = and 
Vc = 0.17, (c) Vs = 0.15 and Vc = 0.08, and (d) V = 0.15 
and Vc — 0.17, respectively. 
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FIG. 4: In-plane optical conductivity as a function of fre- 
quency for different spin (charge) domain-derived scattering 
potentials in the 1/8 antiphase stripe with the SC order pa- 
rameter A = 0. (a) Vs = 0.15 and Vc = 0, (b) T4 = and 
Vc = 0.17, (c) Vs = 0.15 and K = 0.08, and (d) Vs = 0.15 
and Vc = 0.17. 



Rev. Lett. 76, 3412 (1996). 

K. Yamada, C. H. Lee, K. Kurahashi, J. Wada, S. Waki- 
moto, S. Ueki, H. Kimura, Y. Endoh, S. Hosoya, G. Shi- 
rane, R. J. Birgeneau, M. Greven, M. A. Kastner, and Y. 
J. Kim, Phys. Rev. B 57, 6165 (1998). 
B. O. Wells, Y. S. Lee, M. A. Kastner, R. J. Christianson, 
R. J. Birgeneau, K. Yamada, Y. Endoh, and G. Shirane, 
Science 277, 1067 (1997). 

Y. S. Lee, R. J. Birgeneau, M. A. Kastner, Y. Endoh, S. 

Wakimoto, K. Yamada, R. W. Erwin, S.-H. Lee, and G. 

Shirane, Phys. Rev. B 60, 3643 (1999). 

H. A. Mook, P. Dai, S. M. Hayden, G. Aepph, T. G. Per- 

ring, and F. Dogan, Nature 395, 580 (1998). 

T. Valla, A. V. Fedorov, J. Lee, J. C. Davis, and G. D. Gu, 

Science 314, 1914 (2006). 

X. J. Zhou, T. Yoshida, S. A. Kellar, P. V. Bogdanov, E. 
D. Lu, A. Lanzara, M. Nakamura, T. Noda, T. Kakeshita, 
H. Eisaki, S. Uchida, A. Fujimori, Z. Hussain, and Z.-X. 
Shen, Phys. Rev. Lett. 86, 5578 (2001). 



M. Dumm, and D. N. Basov, Seiki Komiya, Yasushi Abe, 
and Yoichi Ando, Phys. Rev. Lett. 88, 147003 (2002). 
M. Ortolani, P. Calvani, S. Lupi, U. Schade, A. Perla, M. 
Fujita, and K. Yamada, Phys. Rev. B 73, 184508 (2006). 
C. Castellani, C. Di Castro, M. Grilh, A. Perah, Physica 
C 341-348, 1739 (2000). 

C. C. Homes, S. V. Dordevic, G. D. Gu, Q. Li, T. Valla, 
and J. M. Tranquada, Phys. Rev. Lett. 96, 257002 (2006). 
Y. Ando, A. N. Lavrov, S. Komiya, K. Segawa, and X. F. 
Sun, Phys. Rev. Lett. 87, 017001 (2001). 
X. J. Zhou, T. Yoshida, A. Lanzara, P. V. Bogdanov, S. A. 
Kellar, K. M. Shen, W. L. Yang, F. Ronning, T. Sasagawa, 
T. Kakeshita, T. Noda, H. Eisaki, S. Uchida, C. T. Lin, 
F. Zhou, J. W. Xiong, W. X. Ti, Z. X. Zhao, A. Fujimori, 
Z. Hussain, and Z.-X. Shen, Nature (London) 423, 398 
(2003). 

M. Dumm, S. Komiya, Y. Ando, and D. N. Basov, Phys. 
Rev. Lett. 91, 077004 (2003). 

M. Sutherland, D. G. Hawthorn, R. W. HiU, F. Ronning, S. 
Wakimoto, H. Zhang, C. Proust, E. Boaknin, C. Lupien, L. 
Taillefer, R. Liang, D. A. Bonn, W. N. Hardy, R. Gagnon, 
N. E. Hussey, T. Kimura, M. Nohara, and H. Takagi, Phys. 
Rev. B 67, 174520 (2003). 

Y. S. Lee, K. Segawa, Y. Ando, and D. N. Basov, Phys. 
Rev. B 70, 014518 (2004). 

T. Tohyama, S. Nagai, Y. Shibata, and S. Maekawa, Phys. 
Rev. Rev. 82, 4910 (1999). 

R. S. Markiewicz, Phys. Rev. B 62, 1252 (2000). 

I. Martin, G. Ortiz, A. V. Balatsky, and A. R. Bishop, 

Europhys. Lett. 56, 849 (2001). 

J. Lorenzana, and G. Seibold, Phys. Rev. Rev. 90, 066404 
(2003). 

A. J. Millis, and M. R. Norman, Phys. Rev. B 76, 
220503(R) (2007). 

P. A. Lee, N. Nagaosa, and X. G. Wen, Rev. Mod. Phys. 
78, 17 (2006). 

J. X. Li, C. Y. Mou, and T. K. Lee, Phys. Rev. B 62, 640 
(2000). 

E. Berg, E. Fradkin, E.-A. Kim, S. A. Kivelson, V. 
Oganesyan, J. M. Tranquada, and S. C. Zhang, Phys. Rev. 
Lett. 99, 127003 (2007). 

The term "insulating" state used here follows Refs. [l3ll6l . 
18 to indicate a strong suppression of Drude peak in the 
optical conductivity. In fact, the spectral weight at the 
Fermi level will not be fully gapped out, so it is not a true 
insulating state. We use the term here is to facilitate our 
comparison with the experiments [Refs. H^ll^lisl ]. 



